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Abstract In the scenario of single-field inflation, this field 
is described in terms of Jacobian elliptic functions. This ap¬ 
proach provides, when constrained to particular cases, ana¬ 
lytic solutions already known in the past, generalizing them 
to a bigger family of analytical solutions. The emergent cos¬ 
mology is analyzed using the Hamilton-Jacobi approach and 
then the main results are contrasted with the recent measure¬ 
ments obtained from the Planck 2015 data. 

PACS 04.20.Fy, 04.20.Jb, 04.40.Nr, 04.70.Bw 


1 Introduction 

One of the main paradigms of the early universe is to the so- 
called inflation, which, for example, allows quantum fluctu¬ 
ations to produce the seeds for the cosmic microwave back¬ 
ground anisotropy and large scale structure, so one of the 
ways by which this stage can emerge is, for example, pro¬ 
tecting the holographic principle [1]. In the standard infla¬ 
tionary cosmology, the universe is dominated by a scalar 
inflaton field, 0, self-interacting through a potential, V{(j)), 
with an energy density = ^0^ -t-V (0), and pressure p,j, = 
i0^ — y (0). Since inflation proceeds if the potential energy 
of the held dominates its kinetic energy, the weak energy 
condition ensures that the pressure is negative during the in¬ 
flationary process allowing for an extremely rapid expan¬ 
sion of the universe. In general, it is not easy to obtain ex¬ 
act solutions of the Einstein equations with an inflaton field 
as source, for a given potential y(0). However, for some 
simple cases, exact solutions can be obtained without even 
making use of the slow-roll approximation. 

To perform a quantitative analysis of the inflation one 
could proceed in two way: using the standard slow-roll ap- 
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proximation [2], or using the Hamilton-Jacobi method [3- 
6]. In this work, we opt for the latter approach, obtaining 
new exact solutions for an inflaton field where the potential 
is expressed in terms of elliptical functions. The reason be¬ 
hind this choice is the fact that many authors have reported 
potentials in terms of trigonometric [7-10] and hyperbolic 
functions[7-15]. Within this context, and keeping in mind 
that both, trigonometric and hyperbolic functions, are noth¬ 
ing but particular cases of elliptical functions, we seek for a 
generalization based on functions of this kind. 

This paper is organized as follows: in Section II we give 
a review of the Hamilton-Jacobi approach to inflation. In 
Section III we present a quartic potential which can be ex¬ 
pressed in terms of elliptical functions, and we provide the 
general solutions and the expressions for the relevant cosmo¬ 
logical parameters. In Section IV we study scalar and tensor 
perturbations and we also compare them with the available 
observational data. 


2 Hamilton-Jacobi approach to inflation 


The starting point is the expression of the equations govern¬ 
ing the background in terms of the inflaton 0, so in units 
where h = c = I these relations are given by 


H^ = 

and 



i02 + y(0) 


( 1 ) 


0-F3H0-Fy'(0)=O, (2) 

where nip = is the Planck mass, a dot represents dif¬ 

ferentiation with respect to the cosmic time, t, and a prime 
represents differentiation with respect to the inflaton scalar 
field. Taking the temporal derivative of Eq. (1), we obtain 

( 3 ) 
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so, using Eq. (2) into Eq. (3) leads to the useful relation 

Also, considering the relation between operators 
d ■ d 

Equation (4) can be written as 

If we assume that the Hubble parameter and the inflaton are 
invertible, then we can write the potential as a function of 
H. Therefore, inserting Eq. ( 6 ) into Eq. (1) we obtain the 
kinematic equation involving the evolution of the scalar field 


(4) 


(5) 


( 6 ) 
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(7) 


c/0 y m'^ 

Since the inflaton is a real scalar field, the condition that 
> pj V must be fulfilled during inflation. Equations 

( 6 ) and (7) are called the Hamilton-Jacobi equations, and 
they dictate the kinematic evolution. In order to obtain ex¬ 
act solutions to Eqs. ( 6 ) and (7), we propose a potential 
that leads to analytic solutions, i.e., that allows us to And 
the generating function, 7/(0). Thus, the set of Eqs. ( 6 ) and 

(7) yield exact solutions. On the other hand, from Eq. ( 6 ), 




a( 0 ) = Ui exp 


a H, it follows that 


An 

in 2 


/ 

J d 


>1’ 7/(0) 


c /0 


(8) 


I, 4 77'(0) 

where a,- = < 2 ( 0 ,), and, from now on, the subscripts i and 
/ imply that quantities are evaluated when inflation begins 
and ends, respectively. The above equation implies that the 
scalar factor a can be written as a function of the inflaton 
field 0 , and thus, assuming that we know the scalar held as a 
function of time, we can obtain the scale factor as a function 
of the cosmological time. 

It is also possible to express the acceleration equation 
for the scale factor as 


-- 772 ( 1 - 6 ;,), 

where the first Hubble hierarchy parameter Eh is given by 
„ 2 \ /^/\2 


Eh = - 


dlnH 


d\na 



(9) 


( 10 ) 


This parameter provides information as regards the acceler¬ 
ation of the universe, so during inflation the bound ^ 1 
is fulfilled, and inflation ends when £„ = 1 . 

Another quantity to describe inflation is the number of 
e-folding of the physical expansion, which is given by 

* ' (11) 


( 12 ) 


A = ln 

so, using Eqs. ( 6 ), ( 8 ) and (10) we obtain 

7'/ (An\ fP' H fP' 1 77' 

N= Hdt = / —-c/0 = /-c/0 

Jti \ ^ A/ 77 


In the description of inflation it is also relevant to show 
that the solutions are independent from their initial condi¬ 
tions. This guarantees the true predictive power that any in¬ 
flationary universe model must have, otherwise the corre¬ 
sponding physical quantities associated with the inflationary 
phase, such as the scalar or tensor spectra, would depend on 
these initial conditions. Thus, any inflationary model needs 
to fulfill the condition that its solutions follow an attractor 
behavior, in the sense that solutions with different initial 
conditions should converge to a unique solution [3]. 

Let us start by considering a linear perturbation, 577(0), 
around a given inflationary solution, expressed by 77(o)(0). 
In the following we will refer to this quantity as the back¬ 
ground solution, and any quantity with the subscript ( 0 ) is 
assumed to be evaluated taking into account the background 
solution. Substituting 77 — 77(o) + 577 into Eq. (7) and then 
linearizing, it is straightforward to show that [ 2 , 16, 17] 


577- 


12 7r 77, 


'PksH'. 


( 0 ) 


which can be solved as 


577(0) = 577(0,) exp 


fP / 3 
\£h 


77' 


77, 


( 0 ) 


c /0 


(13) 


(14) 


where 0 , corresponds to some arbitrary initial value of 0 . 
Since c/0 and 77' have opposite signs (assuming that 0 does 
not change sign due to the perturbation 577) the linear per¬ 
turbations tend to vanish quickly [ 2 ]. 


3 A quartic potential 

In this section we will present a class of potential which 
can be given in terms of elliptic functions. The reason be¬ 
hind this choice is the fact that many authors have reported 
potentials in terms of trigonometric [7-10] and hyperbolic 
functions [7-15]. So, recalling that both functions, trigono¬ 
metric and hyperbolic, are special cases of elliptic functions, 
we will try to And a potential which leads to some function 
of this kind. To proceed, we define the following potential: 
8 ;r 
3 ot2 


y = oo + «! Y{H) + a2Y^{H) -t- a4Y\H), 


(15) 


where the a„, for n= 1,2,3,4, are explicitly given by Oq = 
772 - 772(1 - k2), ai = 277o77i, aj = 2772(1 - k^), a4 = 
H^k^, and the dimensionless function Y is given by 


Y{H) = ^-,^. 
77i 


(16) 


In these definitions, = Hq/Hi, Hi and k are constants 
which will be tested against observational data. 

Substitution of Eq. (15) into the general Eq. (7) implies 

that 




(17) 
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Fig. 1 Plot of the Hubble parameter H as a function of the dimension¬ 
less scalar field <P. 
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Fig. 3 Plot of the parameter of the equation of state (O^ as a function 
of the dimensionless density for different values of the modulus k. 
All plots were made using =1.3. 


where = 12n/m^, and k'^ = 1 — k^. Therefore, perform¬ 
ing an integration (see [18-22]), and then solving for H, we 
obtain 

H{^)=Hi{Ji^+cn[l5{^-^o)]}, ( 18 ) 

where cn(x) = cn(x| k) is the Jacobi elliptic cosine function. 
Fig. 1 depicts //(^) for different values of the modulus k. So, 
by defining the dimensionless scalar field <P ~ (<(> — (po), 

and after a slight manipulation, the potential is conveniently 
written as 

Stt V -IT') 

r(0) = —^— = (J^ + cn<pf-sn^<Pdn^<P, (19) 

3m^ H\ 

where sn(x) = sn(x|A:) is the Jacobi elliptic sine function, 
and dn(x) = dn(x| k) is the Jacobi elliptic delta function. The 
behavior of F(0) is indicated in Fig. 2, where for compar¬ 
ison, different values of the modulus k are shown. The di- 
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Fig. 2 A characteristic behavior of the potential P as a function of the 
dimensionless scalar field <P for different values of the modulus k. All 
plots were made using =1.5. 


mensionless pressure and the energy density are given by 


and 

^ ^ = -(^ + cn4>)Vsn20dn20, (21) 

respectively. Therefore, the parameter of the equation of state, 

(0^ {Pip) = P(p/P(j), can be written as 

(’^(j>{P<j>) = — 1 -f (22) 

where the function fp = f{pp) is given explicitly by 

(l-Jf'2_^2Jf'p'/2_p)[l_;t2(l-jr2-P2Jf'p'/2_p)] 

^2 ■( 23 ) 

Note that the condition f(j, >0 must be satisfied for every 
pair {Jlf,k). 

Also, from Eq. (22) we can see that the Jacobian scalar 
field 0 is driving inflation in such way that the condition 
(0^ < — 1 always is satisfied. Clearly there are four criti¬ 
cal values for the (dimensionless) density which lead to a 
cosmological-like constant regime {(O^ = —1): 

pWp = (J^^±l)2, (24) 

4I4 = (1 - ^) ±2^^^(l-^). (25) 

First of all, p^^ j and p^‘^2 ^re positive and independent of 

the modulus k. But, since 1 > k > 0, then p^ 3 and p^ 4 are 
complex (a conjugate pair) except in the case k =\, where 
they take the value (degenerate) p^^ j = p^‘^2 = see Fig. 

3. 


Using the generating function (18) in Eq. (8), and then 
performing an integration, we find 

at ^ 

where sd(x|k) = sn(x|k)/dn(x|k), and the Gugu function f^(0) 
is given by 


^(0) = cn<(.)_ 


(27) 


sn^* 
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Since this function diverges when k= \ {ox k' = 0) we must 
exclude it and integrate directly this case from Eq. (8), ob¬ 
taining 

fl(0) _ smh0 ( ^(0) , . 


«o 


Qj sinh^>; 
with the function ^(0) defined by 


cosh^> — 1 


= - gcoshcp (29) 

smh5> 

In the same way, using the fact that= —p Hi sn4>dn4> 
together with Eq. (6), we see that the cosmological time is 

f = ln^^"i (0), or more conveniently 

(30) 

if 0 < A: < 1, whereas in the hyperbolic limit. A: = 1, the cos- 

1 

mological time is f = ln,^"i (0), or 

(31) 

Note from Eqs. (26), (30) and (31) we note that a 
gHa(t-u)^ required by the de Sitter stage. 

In addition, making use of the properties of the Jacobian 
elliptic functions, we can write the first Hubble hierarchy 
parameter (10) in the convenient form 
sn^^'dn^0 

[jr + cn0]2- ^ 

In the limit k^Q, this simplifies to 

sin^ ^ 

[jr + COS 0]2’ 

while in the limit k^ 1, we obtain 
sech^^'tanh^^' 

[jr + sech0]2- ^ 

Let us stop briefly in this last expression for the case Jif = 
Hq = 0, which, except for a constant, leads to the expression 
obtained by del Campo [13] when a Chaplygin-like scalar 
field is considered. 

Eor the case of inflation (a/a>Q), the relevant boundary 
condition is de Sitter expansion, £// (0s) = 0. So, for our case 
this condition leads to 0s = 00- From the definition of £//, 
this is equivalent to //'(0) = 0, and it follows that 0o is a 
stationary point of the held; 


= (35) 

From the generating function (18) we see that this point is 
00 = 0. Note that at this point the Hubble parameter is given 
by//(0o)=//o+//i. 

In order to assess whether we have an inflationary pe¬ 
riod, we introduce the deceleration parameter q, defined as 
q = which, in terms of the first Hubble hierarchy pa¬ 
rameter £//, it becomes q = Eh — I, which gives 


q = -l 


sn^ (p dn"^ (P 


(36) 


[,^-|-cn^>]2 

The behavior of q'(0) is indicated in Fig. 4 for different val 
ues of the modulus k. 

The number of the e-folds is found to be 


N{<Pi) = In 


fsd<Pf\ 
\sd0i * 




+ J^\n 


) ■ 


(37) 




Fig. 4 Plot of the deceleration parameter as a function of the dimen¬ 
sionless scalar field <P. 


We must remark that the final value of the dimensionless 
scalar held <Pf is obtained using the transcendental equation 
= —cncl>/± sn^>/^dn^>/^. (38) 

From Eqs. (37) and (38) using N = 60, N = 70, =1.3 

(except for A: = 1 where we have used = 0) and different 
values of the modulus k, some values of are pre¬ 

sented in table 1 . 


The study of the attractor behavior of the model is per¬ 
formed by considering a linear perturbation, 5/7(0), around 
a given solution and following the steps described above. 
Thus, from Eqs. (14), (18), and (32), we get that 


5//(0) = 


sd<Pi /^(0,)\ 

sd0 V^(0) / 




5//(0,). 


(39) 


Since the factor accompanying 5//(0,) is nothing but (a /a;) ^ 
[see the expression (26)], which during inflation increases at 
least 70 e-fold, we can conclude that 5//(0) is very small. 


4 Scalar and tensor pertnrbations 

In the inflationary scenario, the quantum fluctuations are 
found to be relevant because they generate two important 
types of perturbations; density perturbations (arising from 
quantum fluctuations in the scalar field, together with the 
corresponding scalar metric perturbation [23-25]), and relic 
gravitational waves (which are tensor metric fluctuations [26- 
30]). The former is sensitive to gravitational instability and 
leads to structure formation [31], while the latter predicts 
a stochastic background of relic gravitational waves which 
could influence the cosmic microwave background anisotropy 
via the presence of polarization in it [32, 33]. 

In order to describe these perturbations we introduce a 
series of parameters known as the Hubble hierarchy param¬ 
eters. We have already defined one of them through Eq. (10), 
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Table 1 Values for the dimensionless scalar field for some values of the modulus k. The values of were calculated using the condition 
= 1 from Eq. (32), while the values of 4>i were calculated using the condition N{<Pi) = 60 and N{<Pi) = 70 from Eq. (37). In all these 
computations we use = 1.3, except for A: = 1 where we have assumed Jlf = 0. 


1 

1 

1 

A = 60 

1 

1 

A = 70 

1 

k^ 

<Pf 

4>, 



<Pi 


A4> 

0 

1.9515 

1.5764 


0.3752 

1.5757 


0.3760 

0.5 

2.8170 

1.6883 


1.1287 

1.6863 


1.1307 

0.8 

4.0568 

2.0118 


2.045 

2.0088 


2.0480 

0.95 

5.3123 

2.6781 


2.6342 

2.6747 


2.6376 

1 

19.5933 

1.1206 


18.4728 

1.0845 


18.5088 


presented as the first Hubble hierarchy parameter. The sec¬ 
ond Hubble hierarchy parameter, rjH, is defined by 


riH = - 


dlnH' 
d Ina 



(40) 


where plugging in the Jacobian scalar field leads to 
cn0(l-2dn2 0) 


r?// = 


Jf^ + cn^> 


(41) 


The third Hubble hierarchy parameter, is defined by 





(42) 


leading to 

^2 _sn^(Pdn^(I>{6k^sn^^-4k--l) 

[jr+cn0]2 ■ ^ 

The quantum fluctuations produce a power spectrum of 
scalar density fluctuations of the form [34, 35] 


-V 

27tJ 


(44) 


aH=k 


This perturbation is evaluated when a given mode k crosses 
outside the horizon during inflation, i.e. at aH — k. These 
modes do not evolve outside the horizon, so we can assume 
they keep a fixed value after crossing the horizon during in¬ 
flation. In order to obtain some comparison with the avail¬ 
able observational data, we introduce the scalar spectral in¬ 
dex Hs defined as 
d\n^^ 

ns-l = -(45) 

c/lnk 

After a brief calculation, one obtains that 
ns = 1— 4e// -\-2rjH, (46) 

so, using Eqs. (32) and (41), the spectral index becomes 
2|,jf’cn4>(i-2dn^ <!>)+I^cn^ <f>-3dn^ if>(|cn-<h-i)j I 
\j¥+cn^ • 

It is well known, that not only scalar curvature pertur¬ 
bations are generated during inflation. In addition quantum 
fluctuations generate transverse-traceless tensor perturbations 
[31], which do not couple to matter. Therefore they are only 
determined by the dynamics of the background metric. The 
two independent polarizations evolve like minimally cou¬ 
pled massless fields with spectrum 
1 ^ ^ ^ 
m2. 


-V 

27tJ 


(48) 


aH=k 



Primordial tilt (rts) 

Fig. 5 This plot shows the parameter r as a function of the scalar spec¬ 
tral index for various values of the modulus k, i.e. k = 0,ld = 0.5, 
id = 0.8, k^ = 0.95 and k — 1. The marginalized joint 68% and 95% 
confidence level regions using Planck TT+ lowP or Planck TT, TE, EE 
+ lowP are shown [36]. Also, constraints from the Planck 2013 data 
release are also shown for comparison. Every curves have been plotted 
with H\ = I, = 1.3, and they are normalized in such way that r = 0 
when Its = i. 


In the same way as the scalar perturbations, it is possible to 
introduce the gravitational wave spectral index iij defined 
by tiT = ^ which in our case becomes nj = —2e// . At 

dink 

this point, we can introduce the tensor-to-scalar amplitude 
ratio r = which becomes 

r = AeH, (49) 

and, in our case, we obtain that 
4sn^4>dn^^' 

[jr + cn0]2- 


In Fig. 5 we display the parameter space of r, given by Eq. 
(50), and the spectral index given by Eq. (47). 


5 Final remarks 

There is no doubt that elliptic functions are a powerful tool 
to describe natural phenomena. In this paper we have used 
them to describe a scalar field in order to find inflation¬ 
ary solutions. In particular, we have used an inflationary 
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universe model in which the inflaton field is characterized 
by a parameter of the equation of state co^ = — 1 + -s/fy, 
where fy is given by Eq. (23), and its kinematical evolution 
was described by the Hubble parameter given by //(0) = 
+ cn [j3 (0 — 0o)]}. A detailed analysis shows that 
the allowed values for //(0) depend strongly on the values 
of the pair On this basis, the scalar potential, V (0), 

the corresponding number of e-foldings, and the attractor 
feature of the model were described. 

We should mention here that the first hierarchy parame¬ 
ter £// (and thereby, the deceleration parameter q) related to 
the inflaton field results in such a way that it is possible to 
reproduce a generalized Chaplygin gas by making 0. 

This analogy encouraged us to start a generalization of the 
GCG in terms of elliptic functions [37]. 

It was also shown that the potential works quite well 
when compared against the measurements recently released 
based on Planck data. This situation is the main motivation 
to study inflationary universe models with this kind of scalar 
potential. 

Finally, in general terms, we have found that the tensor- 
to-scalar ratio can adequately accommodate the currently 
available observational data for some values of the 
parameters. In this context, we have shown that the model 
here presented is appropriate for describing inflationary uni¬ 
verse models. 


Appendix A: A brief review of Jacobian elUptic 
functions 


As a starting point, let us consider the elliptic integral 
u(y,k) = u = f ■ 

Jo 


r<P 

Jo 


=F{(p,k), 


, , -—^ - vr,-., (A.l) 

-'0 vl — sin^ 9 
where F{(p,k) is the normal elliptic integral of the first kind, 
and k is the modulus. The problem of the inversion of this 
integral was studied and solved by Abel and Jacobi, and this 
leads to the inverse function defined by y = sin^ = sn(M,k) 
with (p — ami/, and the functions are called Jacobi elliptic 
sine u and amplitude u. 

The function sn u is an odd elliptic function of order 
two. It possesses a simple pole of residue 1 /k at every point 
congruent to iK' (mod 4K, 2iK') and a simple pole of residue 
— 1 /k at points congruent to 2K -|- iK' (mod 4 K, 2iK'), where 
K = K{k) = F{7t/2,k) is the complete elliptic integral of the 
first kind, K' = F{K/2,k'), and k' = s/l—kf is the comple¬ 
mentary modulus. 

Two other functions can then be defined by cn(M,k) = 
\/l — y2 = cos 9 , dn(i/,k) = a/I — k^y^ = A(p = 
—k^ sin^. The set of functions { sn i/, cn m , dn i/} are called 


Jacobian elliptic functions, and they take the following spe¬ 


cial values: 

sn(i/|0) = sini/, sn(i/|l) = tanhi/, (A.2) 

cn(M|0) = cost/, cn(i/|l) = sechi/, (A.3) 

dn(i/|0) = l, dn(i/|l) = sechi/. (A.4) 

Some fundamental relations between Jacobian elliptic func¬ 
tions are 

sn^M-f cn^M = 1, (A.5) 

k^sn^i/-f dn^M = 1, (A.6) 

dn^M — k^cn^M = k'^, (A.7) 

k'^svJu -f cn^M = dn^M. (A.8) 
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